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ABSTRACT 


A relationship between a person's lean body mass and the amount of maximum 
torque that can be produced with each isolated joint of the upper extremity was 
investigated. The maximum dynamic isolated joint torque (upper extremity) on 14 
subjects was collected using a dynamometer multi-joint testing unit. These data were 
reduced to a table of coefficients of second degree polynomials, computed using a 
least squares regression method. All the coefficients were then organized into look-up 
tables, a compact and convenient storage/retrieval mechanism for the data set. Data 
from each joint, direction and velocity, were normalized with respect to that joint's 
average and merged into files (one for each curve for a particular joint). Regression 
was performed on each one of these files to derive a table of normalized population 
curve coefficients for each joint axis direction and velocity. In addition, a regression 
table which included all upper extremity joints was built which related average torque 
to lean body mass for an individual. These two tables are the basis of the regression 
model which allows the prediction o^ dynamic isolated joint torques from an 
individual's lean body mass. 
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1.0 INTRODUCTION 


A relationship exists between a person's lean body mass and the amount of maximum 
torque that can be produced with each isolated joint of the upper extremity. The use of 
an easily measured parameter (lean body mass) to predict dynamic isolated joint 
torque, which is time consuming to measure, would be extremely valuable. In this 
study, we have collected maximum dynamic isolated joint data on 14 subjects for the 
upper extremity and formulated a regression model which will allow prediction of a 
joint angle versus torque curve for a particular individual. There are three phases in 
our study: 1) data collection, 2) data reduction, and 3) model formulation. 

Phase one of our study, data collection, involved measuring the maximum torque for 
all the upper extremity joints (shoulder, elbow, and wrist) at 4 velocities for 1 4 subjects. 
Data for three subjects' lower extremities (hip, knee, and ankle) were also collected. 
Torque was measured by using a LIDO multi-joint testing unit (Loredan Biomedical, 
Inc., West Sacramento, California). The subjects were positioned so that the axis of the 
joint was directly in line with the axis of the dynamometer goniometer. Dynamometer 
attachments were selected and used to isolate the joint being measured. In this 
manner, all joints were characterized for all axes of rotation over a range of velocities. 
The data were collected using the Loredan software, LIDO Active 3.3, on an IBM PC. 
For all cases, the data set consisted of torque and angle pairs. In addition, 
anthropometric data were also collected which included height, weight, age, sex, skin 
fold measures, and dimensional assessment according to the format specified in 
NASA Man-Systems Integration Standards (MSIS) document [1]. 

The second phase of the project, data reduction, began with transferring the data to a 
UNIX-based workstation (Silicon Graphics). These data were formatted into an ASCII 
file, noise filtered, reformatted, and reduced to a table of coefficients of second degree 
polynomials. The polynomial coefficients were computed using a least squares 
regression method. These polynomials represent the torque as a function of angle 
(i.e., torque = a+ b*angle + c* angle**2, where a, b, and c are the polynomial 
coefficients). All the coefficients were then organized into look-up tables. These 
tables represent a compact and convenient storage/retrieval mechanism for our entire 
data set and are available upon request. 

The third phase of our project involved model formulation. Data from each joint, 
direction and velocity, of an individual were normalized with respect to that joint's 
average and merged into files (one for each curve for a particular joint). Regression 
was performed on each one of these files to derive a table of normalized population 
curve coefficients for each joint axis, direction, and velocity. In addition, a regression 
table, which included all upper extremity joints, was built which related average torque 
to lean body mass for an individual. These two tables are the basis of the regression 
model which allows isolated joint curve prediction. Because of the limited number of 
subjects (3) for the lower extremity data set, no correlation to lean body mass is 
presented here for these measurements. 
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Finally, we have encapsulated the results of this study in a tool kit of software routines 
executing on a variety of platforms such as UNIX, DOS, and Macintosh machines 
(Appendix D). This code contains our isolated joint torque model (with all the tables of 
torque coefficients) and will allow the prediction of dynamic isolated joint torques from 
an individual's lean body mass. 
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2.0 OBJECTIVE 


Prediction equations have been developed from lean body mass for isometric/static 
strength [2]. Few studies have utilized lean body mass to predict isokinetic or dynamic 
strength [3]. An absence of literature exists when correlating lean body mass to 
isolated joint isokinetic mean torque or torque over an entire range of joint motion. In 
addition, lean body mass has not been used to predict the mean torque for dynamic 
complex tasks involving multiple joints. The use of an easily measured parameter 
(lean body mass) to predict dynamic isolated joint torque would be of significant 
value. 

The objective of this project is to develop prediction equations which can be used to 
calculate isolated joint torque (either mean torque or torque as a function of angle) and 
mean torque for a complex task from the measurement of a person's lean body mass. 

Specific aims: 

1 . Document all data measurement and data processing techniques. 

2. Develop prediction equations for: 

a. Mean dynamic torque over an entire range of motion for a particular 
joint (shoulder, elbow, and wrist) and during a complex task (ratchet 
wrenching or an extravehicular task). 

b. Dynamic torque-position curves over the entire joint range of motion 
from lean body mass for individual isolated joints (shoulder, elbow, and 
wrist). 

3. Develop a set of software tools which simplify dynamic torque data access, 
manipulation, and prediction for the set of isolated joints measured. 
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3.0 METHOD 


3.1 Data Collection 
Subjects 

Fourteen subjects (8 males and 6 females) aged 21 to 28 years volunteered and 
participated in this investigation. The study was evaluated and approved by the 
institutional Review Board of the University of Texas Medical Branch at Galveston, 
Texas. All subjects were informed of potential risks and signed a consent to 
participate in the study. 

Equipment 

The LIDO multi-joint testing unit (Loredan Biomedical, Inc., West Sacramento, 
California, Figure 1) is an integrated system consisting of a dynamometer connected to 
a personal computer via RS232c lines. The force unit comes with a series of 
attachments and a subject bench which allows isolation of particular joints. The 
software allows precise control of the actuator head and the various modes of 
operation (e.g., isometric, isokinetic, and concentric). In addition, a database of 
subjects can be maintained and displayed with the provided graphical software. The 
unique feature of this system is that it outputs all data in a machine-readable form for 
more accurate data analyses. This system was used to measure all the isolated joint 
forces as well as the composite test cases. 

Data Collection Procedures 

Data collection for this project occurred over an 8-week period at the University of 
Texas Medical Branch, School of Allied Health Sciences, Human Performance 
Laboratory. The general procedure for evaluating all the upper and lower extremity 
joint movements and simulated ratchet maneuvers was the same. Specific subject 
and joint positioning for the isometric and isokinetic tests is described in each 
independent section. The time required to perform all of the dynamic and isometric 
shoulder, elbow, and wrist joint measurements was 1 .5 to 2 hours; evaluation of the 
hip, knee, and ankle measures required 1 .0 to 1 .5 hours; and evaluation of the 
3 ratchet wrench maneuvers required 0.5 to 1 .0 hours per subject. 
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Actuator 



Controller 


Figure 1. LIDO multi-joint testing unit. 

On each testing day subjects reported to the laboratory in a fasting condition (food was 
restricted 3 hours prior to evaluation). For all upper extremity tests, the subject was 
positioned so that the axis of the joint was directly in line with the axis of the 
dynamometer goniometer. All measures were taken without gravity compensation. 
Maximum isometric contraction (MIC) measures were taken first. Dynamometer 
attachments were selected and placed to isolate the joint; the subject was positioned 
on the instrument and maximally stabilized; and then the joint was positioned at a 
specific angle. The subject was instructed each time to give maximum efforts. The 
subject performed three submaximal contractions followed by one MIC. The subject 
was given 3 minutes of rest and this procedure was repeated in the opposing 
direction. Next, the subject recovered for 5 minutes and then performed the isokinetic 
testing. The subject was instructed to give maximum efforts for each repetition and to 
move through the entire range of motion as rapidly and as forcefully as possible. The 
4 joint velocity settings (60, 120, 180, and 240 deg/sec) were randomly assigned. The 
subject performed three submaximal contractions at the designated velocity followed 
by five maximum contractions. The subject was given 3 minutes of recovery between 
each velocity setting. 


Shoulder Flexion and Extension 

Subject was in a supinated position and was stabilized with velcro straps at the waist 
and across the clavicle just proximal to the shoulder to keep the right shoulder firmly in 
contact with the plinth (Figure 2). The angle for isometric testing was 90 degrees. A cuff 
attached to the dynamometer arm was placed on the upper arm just proximal to the 
elbow joint. Therefore, the shoulder joint was isolated and force was applied at the 
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point of the cuff attachment. The range of motion where shoulder flexion and extension 
torque was measured was between 20 and 180 degrees of shoulder flexion. 

Shoulder Abduction and Adduction 

Subject was lying on the side and was stabilized with velcro straps at the pelvis and 
across the upper chest at the axilla line to keep the subject's chest wall firmly in 
contact with the plinth (Figure 3). The angle for isometric testing was 90 degrees. A 
cuff attached to the dynamometer arm was placed on the upper arm just proximal to 
the elbow joint. Therefore, the shoulder joint was isolated and force was applied at the 
point of the cuff attachment. The range of motion where shoulder abduction and 
adduction torque was measured was between 15 and 145 degrees of shoulder 
abduction. 



Figure 2. Shoulder flexion and extension. 
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Figure 3. Shoulder abduction and adduction. 


Shoulder Internal (Medial) and External (Lateral) Rotation 

Subject was in a supinated position with the shoulder placed at 90 degrees of 
abduction. The subject was stabilized with velcro straps at the waist and across the 
mid-upper arm to keep the right-upper arm and elbow firmly in contact with the plinth 
(Figure 4). The angle for isometric testing was 0 degrees. A cuff attached to the 
dynamometer arm was placed on the forearm just distal to the elbow joint. The 
shoulder joint was isolated and force was applied at the point of the cuff attachment. 
The range of motion where shoulder internal and external rotation torque was 
measured was between 0 and 70 degrees of internal rotation and 0 and 60 degrees of 
external rotation. 

Elbow Flexion and Extension 

Subject was in a supinated position and was stabilized with velcro straps at the waist 
and across the chest and mid-upper arm to keep the right-upper arm and elbow firmly 
in contact with the plinth (Figure 5). The angle for isometric testing was 60 degrees. A 
cuff attached to the dynamometer arm was placed on the forearm just proximal to the 
wrist joint. The elbow joint was isolated and force was applied at the point of the cuff 
attachment. The range of motion where elbow flexion and extension torque was 
measured was between 15 and 135 degrees of elbow flexion. 
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Wrist Flexion and Extension 

Subject was in a sitting position and was stabilized with velcro straps at the waist and 
across the chest to keep the subject's back firmly in contact with the seat (Figure 6). 

An attachment for the forearm was secured to the side of the LIDO. The subject's right 
forearm was secured in a supinated position into the device by three small velcro 
straps: 1) just distal to the elbow joint, 2) mid-forearm, and 3) just proximal to the wrist 
joint (Figure 7). This configuration was designed to keep the forearm firmly in contact 
with the forearm stabilizing device. The angle for isometric testing was 0 degrees. 

The subject gripped a handle device attached to the dynamometer shaft. The wrist 
joint was isolated and force was applied at the point of the handle attachment. The 
range of motion where wrist flexion and extension torque was measured was between 
0 and 60 degrees for wrist flexion and 0 and 45 degrees for wrist extension. 



Figure 4. Shoulder internal (medial) and external (lateral rotation). 
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Figure 6. Wrist flexion and extension. 


9 


ORIGINAL PAGE 

BLACK AND WHITE. PHOTOGRAPH 


Figure 7. Closeup view of wrist flexion and extension. 


Wrist Radial and Ulnar Deviation 

Subject was in a sitting position and was stabilized with velcro straps at the waist and 
across the chest to keep the subject's back firmly in contact with the seat (Figure 6). 

An attachment for the forearm was secured to the side of the LIDO. The subject's right 
forearm was secured in a supinated position into the device by three small velcro 
straps: 1) just distal to the elbow joint, 2) mid-forearm, and 3) just proximal to the wrist 
joint (Figure 8). This configuration was designed to keep the forearm firmly in contact 
with the forearm stabilizing device. The angle for isometric testing was 0 degrees 
(neutral). The subject gripped the handle device attached to the dynamometer shaft. 
Therefore, the wrist joint was isolated and force was applied at the point of the handle 
attachment. The range of motion where wrist radial and ulnar deviation torque was 
measured was between 0 and 35 degrees for wrist radial deviation and 0 and 
35 degrees for wrist ulnar deviation. 

Wrist (Forearm) Pronation and Supination 

Subject was in a sitting position and was stabilized with velcro straps at the waist and 
across the chest to keep the subject's back firmly in contact with the seat (Figure 6). 

An attachment for the forearm was secured to the side of the LIDO. The subject's right 
forearm waa in a neutral position and 2 small velcro straps were placed: 1) just distal 
to the elbow joint and 2) mid-forearm (Figure 9). This configuration was designed to 
keep the proximal section of the forearm firmly in contact with the forearm stabilizing 
device and to allow pronation and supination to occur. The angle for isometric testing 
was 0 degrees (neutral). The subject gripped the handle device attached to the 
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dynamometer shaft. Therefore, the forearm was isolated and force was applied at the 
point of the handle attachment. The range of motion where wrist pronation and 
supination torque was measured was between 0 and 60 degrees for wrist pronation 
and from 0 to 60 degrees for wrist supination. 



Figure 8. Wrist radial and ulnar deviation. 
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Figure 9. Wrist (forearm) pronation and supination. 


Hip Flexion and Extension 

Subject was in a supinated position and was stabilized with velcro straps at the chest 
and waist and across the left-upper thigh to keep the back and pelvis firmly in contact 
with the plinth (Figure 10). The angle for isometric testing was 90 degrees. A large 
cuff attached to the dynamometer arm was placed on the thigh just proximal to the 
knee joint. Therefore, the hip joint was isolated and force was applied at the point of 
the cuff attachment. The range of motion where hip flexion and extension torque was 
measured was between 0 and 1 1 0 degrees of hip flexion. 

Hip Abduction and Adduction 

Subject was lying on the side and was stabilized with velcro straps at the chest and 
waist and across the left-upper thigh to keep the subject's chest wall and pelvis firmly 
in contact with the plinth (Figure 11). The angle for isometric testing was 0 degrees. A 
large cuff attached to the dynamometer arm was placed on the thigh just proximal to 
the knee joint. Therefore, the hip joint was isolated and force was applied at the point 
of the cuff attachment. The range of motion where hip abduction and adduction torque 
was measured was between 0 and 60 degrees of hip abduction. 

Hip Internal (Medial) and External (Lateral) Rotation 

Subject sat on an elevated table and was stabilized with velcro straps at the waist and 
across the left-upper thigh to keep the subject's pelvis firmly in contact with the 
exterior table. The angle for isometric testing was 0 degrees (neutral). A cuff attached 
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to the dynamometer arm was placed on the lower leg just distal to the knee joint. 
Therefore, the hip joint was isolated and force was applied at the point of the cutt 
attachment. The range of motion where hip internal and external rotation torque was 
measured was between 0 and 15 degrees of internal rotation and from 0 to 
35 degrees of external rotation. 



Figure 10. Hip flexion and extension. 



Figure 11. Hip abduction and adduction. 
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Knee Flexion and Extension 


?^® Ct ^ as J n ® sit j ln 9 position and was stabilized with velcro straps at the waist and 
acrose the chest to keep the subject's back firmly in contact with the seat (Fiqure 12). 

i inrf U Tu- CtS n ? bt thl 9 h was stabi | iz ®d with an attachment secured to the side of the 
stahnwln Hm?i^ 9Ur ?I£ n Was . designed to keep the thigh firmly secured between the 
atShoH? d fh f and the Seat ' The ang,e for isometric testing was 60 degrees. A cuff 
an^io Th d y namo "^ er arm , was placed on the lower leg just proximal to the 

_ f +ha r,iff oJi!?h ref0r ?’ t? knee J0ln ] was isolate d and force was applied at the point 
attachme ^t- The range of motion where knee flexion and extension torque 
was measured was between 5 and 100 degrees of knee flexion. 

Ankle Plantarflexion and Dorsiflexion 


Sokx 01 3 , su P |nated Position with the knee joint fully extended. The subject was 
stabilized with velcro straps at the chest and waist to keep the subject's back and 
peivis firmly in contact with the plinth (Figure 13). The subject's right knee was 

* W, l 1 an fttachnient secure d to the side of the LIDO. This configuration was 
2 T d ‘° ke ®P the knee joint firmly secured between the stabilizing device and the 

Slc th ^ The ^ a H 9 ?u for lsomatric testing was 0 degrees (neutral). The subject's right foot 
was secured by three small velcro straps and placed into a foot plate device 9 

thi 9 olw • ■ T he P late device was attached to the dynamometer shaft. Therefore, 
the ankle joint was isolated and force was applied at the foot plate attachment. The 

ran ge of motion where ankle plantarflexion and dorsiflexion torque was measured 

W3 m batwe ® n °. and 30 de grees for ankle plantarflexion and from 0 to 10 degrees for 
ankle dorsiflexion. a 



Figure 12. Knee flexion and extension. 
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Figure 13. Ankle plantarflexion and dorsiflexion. 




Figure 14. Closeup view of ankle plantarflexion and dorsiflexion. 
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Simulated Ratchet Pushing and Pulling Maneuver Torque 

The axis of the dynamometer goniometer was positioned in line with the subject's 
greater trochanter. For all simulated ratchet pushing and pulling maneuvers, the 
subject was in a sitting position. The right upper extremity was evaluated. The subject 
was stabilized with velcro straps at the waist and across the chest to keep the subject's 
back firmly in contact with the seat (Figure 15). The subject gripped a simulated 
ratchet device at a height of 90% of the linear distance measured from the subject's 
greater trochanter to the acromioclavicular joint (Figure 16). The subject was 
instructed to move the ratchet device from a full back position (shoulder joint at 
10.9 +/- 0.8 degrees; elbow joint at 1 12.0 +/- 6.0 degrees) to a full forward position 
(shoulder joint at 66.0 +/- 6.8 degrees; elbow joint at 0 degrees) [ratchet push 
maneuver], and then to return to the starting full back position [ratchet pull maneuver]. 
The subject was told to give maximum efforts for each repetition and to move through 
the entire range as rapidly and forcefully as possible, but to not allow the scapula 
(upper back) to lose contact with the seat. Individual subject measures of ratchet grip 
height and shoulder and elbow joint angle of excursion were recorded. For all 
ratcheting conditions, the subject performed three submaximal contractions followed 
by five maximum contractions. The subject was given 5 minutes of recovery between 
each of the 3 ratcheting conditions. The first ratcheting condition consisted of two 
trials. The first trial consisted of a constant resistance of 25 ft-lbs (isotonic load) for the 
ratchet push maneuver and 0 ft-lbs of resistance for the pull manuever. The subject 
was given a 5-minute recovery period and then performed the second trial with 0 ft-lbs 
of resistance for the ratchet push maneuver and a constant resistance of 25 ft-lbs for 
the pull maneuver. The second condition was the push and pull ratcheting done at a 
constant angular velocity of 120 deg/sec (isokinetic). The third condition was the push 
and pull ratcheting done at a constant angular velocity of 240 deg/sec (isokinetic). For 
each repetition, torque and angle position data were measured. 
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Figure 15. Front view of simulated ratchet pushing and pulling maneuver. 



Figure 16. Side view of simulated ratchet pushing and pulling maneuver. 
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Assessment of Anthropometric Data 


Anthropometric data were gathered for height, weight, age, sex, skin fold measures, 
and dimensional assessment. Height was measured in centimeters and weight in 
kilograms on a physician scale. Skin fold measures were taken with a Lange skin fold 
caliper. Males were measured at the abdomen and anterior thigh and chest; and 
females were measured at the anterior thigh, supraillac crest, and triceps. Estimation 
of body fat and lean body mass were performed using equations specific for sex and 
age [4]. See Table 1 for a summary of this data. 


TABLE 1 . Lean Body Weight 


SUBJ 

SEX 

AGE 

HT (cm) 

WT (kg) 

%Body Fat 

Lean Body 

- 

1. 

M 

23 

173 

64.9 

10.4 

58.2 


2. 

M 

25 

178 

76.0 

10.4 

68.1 


3. 

M 

28 

188.3 

84.5 

14.5 

72.2 


4. 

F 

23 

172.3 

7.2 

29.5 

54.4 

- 

5. 

M 

22 

185.0 

88.2 

10.7 

78.8 

- 

6. 

M 

25 

180.0 

86.0 

5.7 

81.1 


7. 

M 

26 

176.3 

95.2 

19.0 

77.1 

- 

8. 

F 

23 

174.0 

60.6 

18.6 

49.3 

- 

9. 

F 

22 

168.0 

59.1 

14.8 

50.4 


10. 

F 

21 

158.5 

46.4 

17.2 

38.4 

1 

11. 

F 

21 

158.0 

51.6 

19.5 

41.5 

s 

12. 

M 

21 

178.5 

80.7 

8.9 

73.5 

- 

13. 

F 

23 

166.0 

55.7 

20.8 

44.1 

= 

14. 

M 

23 

162.5 

68.5 

11.3 

60.8 



The dimensional anthropometric data were taken with cloth tape measures [1]. The 
time required to make all anthropometric measures was 0.25 to 0.5 hours per subject. 


3.2 Data Reduction 

A set of streamlined programs was developed to process the raw strength data 
(produced directly by the LIDO force torque dynamometer) into a compact polynomial 
coefficient format. The raw data were collected using the LIDOACT software executing 
on an IBM PC. The files produced on the PC were transferred to the VAX system using 
a data communication software package (Kermit) in binary mode. These data were 
then transferred to a UNIX-based workstation. 

After being separated into files by subject, velocity, direction, and degree of freedom 
for each joint, each torque versus angle data file was viewed graphically and edited for 
extraneous data points (Figure 17). Figure 17 shows that the initial and final portions 
of the curves were omitted because of the startup time during which the subject is 
beginning to apply a maximum torque. At the end of a motion, the subject was 
anticipating the stopping and change of direction of the LIDO actuator arm. These 
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transition regions of torque were inconsistent and so were not part of our modeling 
effort. 


Elbow Flexion 


Raw Data 



Filtered Data 

40.0| 



o.o I . . . . . 1 

0.0 75.0 150.0 

Degrees 


Figure 17. Data before and after visual editing. 


The following is the flow of data from raw files to a torque function coefficient file 
(Figure 18). This is the actual UNIX script file used: 

#uncompress the large data files 
uncompress *.raw 

"Is" sh*.raw| lido | tosort | sort | toffc > right_shoulder.ffc 
"Is" el*.raw| lido | tosort | sort | toffc > right_elbow.ffc 
“Is” wr*.raw| lido | tosort | sort | toffc > right_wrist.ffc 
#clean up 
compress *.raw 
mv *.xy* xy 
rm *.asc 

The following is an explanation of the above script file: 

1 . The first line of the script file uncompresses the data files needed. 

2. The "Is" command feeds a list of the file names (one at a time) to the LIDO 
program. 
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3. The LIDO program converts each file it receives from the LIDO format into an 
ASCII format and passes the data (one file at a time) to tosort . 

4. The tosort program computes the regression equations and collates person 
data and passes that condensed data on to the sort program. In addition, it 
creates xy files of force versus angle and stores them for review later. 

5. The sort program sorts the data on each field and passes that on to the toffc 
program. 

6. The toffc program processes sorted data to produce files in the torque 
function coefficient format (Figure 18). 

7. The original data is then recompressed in the seventh line to conserve disk 
space. 

8. The xy files are moved into a separate directory called xy in the eighth line. 

9. All the unnecessary files are cleaned up in the ninth line. 

The following is a quick description of each of the processing programs: 

lido.c - This program converted a raw LIDO generated data file to a standard 
output ASCII file. The LIDO generated data file (extension .raw) had a main 
header area and several subheader areas. Lines in these areas were terminated 
with a line feed as in a normal ASCII text file. The data lines were terminated with 
a carriage return and no line feed. The data lines were three integer values of 
four characters each. The first value was the corrected torque at the cuff [5, 6], the 
second value was the angle, and the third was the uncorrected torque at the shaft 
of the LIDO actuator. The corrected torque values were chosen for further 
processing. 

tosort.c - This program converted the ASCII LIDO data file to torque function 
coefficients and an associated subject name, joint name, direction, and velocity. 
The output was used as standard input for the UNIX sort utility. This sorted output 
was then piped to the toffc.c program. 

toffc.c - This was a program which processed output from the tosort program and 
converted it into a torque coefficient file (Figure 18). 
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in nun mm im mm in inn him mu n i 


right_shoulder /* joint name*/ 
x /* axis */ 

abduction /* direction*/ 

4 /*number of velocities*/ 

/* ^ velocity , 3 polynomial coefficients, note: Y = A + Bx + Cx**2 where A.B.C are the coefficients*/ 


60.000000 

120.000000 

180.000000 

240.000000 
adduction 

4 

60.000000 

120.000000 
180.000000 

240.000000 

y 

extension 

4 

60.000000 

120.000000 
180.000000 

240.000000 
flexion 

4 

60.000000 

120.000000 
180.000000 

240.000000 
z 

lateral 

4 

60.000000 

120.000000 
180.000000 

240.000000 
medial 

4 

60.000000 

120.000000 
180.000000 
240.000000 


3.477892E+01 
5.058879E+01 
3.441 185E+01 
4.630580E+01 


3.899835E+01 
3.210880E+01 
1.1 2871 2E+01 
2.879990E+01 


1 .853307E+01 
1.598296E+01 
1.51 601 0E+01 
9.10331 1E+00 


6.689232E+01 

4.769003E+01 

5.670030E+01 

4.729467E+01 


1 .650289E+01 
1 .434598E+01 
1 .8731 10E+01 
1 .489092E+01 


2.395702E+01 
2.124104E+01 
2.144591 E+01 
1 .980874E+01 


2.220639E-01 

-1.974063E-01 

-2.165272E-02 

-2.306546E-01 


-2.004480E-02 
6.244854E-02 
4.01 3751 E-01 
-4.25531 2E-02 


1 .233500E-01 
2.713108E-01 
2.109192E-01 
1. 32561 5E-01 


-4.941 486E-01 
-2.906334E-01 
-4.014365E-01 
-3.183639E-01 


-4.517265E-02 
-2.340409E-01 
-1 .003202E-01 
-4.907729E-02 


1.128039E-01 
1 .564447E-01 
1. 82661 7E-01 
1. 33621 7E-01 


-2.821 324E-03 
-2.446309E-04 
-1.170670E-03 
6.999267E-05 


-1.001798E-04 
-3.269533E-04 
-1 .763708E-03 
6.87231 5E-04 


-1.150867E-04 
-9.97541 8E-04 
-8.056303E-04 
-4.60291 2E-05 


1 .258664E-03 
8.198689E-04 
1.149240E-03 
8.940246E-04 


-5.588200E-04 
-2.806892E-03 
-1 .389751 E-03 
4.31 1464E-04 


-1 .947907E-04 
-3.186225E-03 
-1 .797236E-03 
-8.705539E-04 


Figure 18. Example of the coefficient file format. 
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3.3 Begression Model De velopment 

Creating the regression model involved processing the data from each joint, direction 
and velocity, into separate files. There were three phases in the processing of these 
files which allow the prediction of individual torque curves. 

1 . Creating lean body mass to average torque relationship tables. 

2. Generating a population shape coefficient table. 

3. Linking the population shape table to the lean body mass average torque 
table. 

Lean Body Mass to Average Torque Table Generation 

A relationship exists between lean body mass and average torque. Figure 19 shows 
representative examples of this linear relationship for the shoulder, wrist, and elbow. 
Correlation coefficient values range from 0.80 to 0.95. The lean body mass for each 
of the subjects was calculated from the measured weight and % body fat: 

lean body mass = weight * (1 - % body fat) 

The average torque was then computed for the entire joint range for each velocity of 
each joint. From these data, linear regression coefficients were obtained and 
organized into coefficient tables. These tables can be used to calculate the average 
torque for a particular joint and velocity when the lean body mass is known. 

Obtaining Population Curves 

In our model an assumption was made that for an isolated joint motion, individuals use 
the same muscle groups. These muscle groups basically have the same shape, 
orientation, and points of attachment and differ in the magnitude of force exerted. 
Hence, when the torque-position curves are normalized and plotted, a general trend 
can be seen for a particular joint. Figure 20 represents the normalized data for all the 
subjects for the elbow and shoulder flexion. These data show a definite trend. To 
take advantage of this relationship, the data were processed in the following way: 

1 . We normalized each individual data set (isolated joint for each axis, 
direction, and velocity) with respect to the average torque for that data set. 

2. For all subjects, the normalized data for an isolated joint (axis, direction, 
and velocity) were combined into one file. A total of 7 axes, 2 directions, 
and 4 velocities resulted in 56 composite files of normalized subject data. 
(Figure 20 shows a sample case.) 

3. A regression for each of these combined files was then computed. 

4. The second order polynomials, representing the shapes of the joint motions, 
were then organized into coefficient tables (Figure 18). These tables allow 
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the computation of a normalized torque-position curve for any joint, axis, 
direction, and velocity. 

Generation of an Individual's Torque-Position Curve 

To generate a torque-position curve from lean body mass information it is necessary to 
multiply the normalized torque-position curve coefficients (established in the table from 
step 2) by the average torque. The average torque is computed from the data in the 
lean-body-mass to average-torque-coefficient table established in step 1 above. 

These two tables provide the data to compute a polynomial which represents an 
individual's torque-position curve for that joint, axis, direction, and velocity. 

For all of our data, torque was measured for four separate velocities. No attempt was 
made to find a correlation between these velocities. Velocity dependence was not 
consistent among all joints. Since dynamic strength is important in most joint motions, 
we feel that this area requires more attention. 
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Shoulder Flexion 


Shoulder Abduction 



v(1-%bf) v(1-*bf) 


Shoulder Lateral Rotation Elbow Flexion 



30.0 60.0 90.0 30.0 60.0 90.0 

wO-'SCbf) v(1-*bf) 


Figure 19. Lean body mass (kg) versus mean torque (ft-lbs). 
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Wrist Flexion 


Wrist Radial Deviation 



wO-'SEbf) v(1-*bf) 


Wrist Supination 



v(l-%bf) 


Figure 19. Lean body mass (kg) versus mean torque (ft-lbs) (continued). 
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Elbow Flexion: velocity = 120 deg/sec 
Normalized Date: 1 4 Subjects Polynomial Fit to Normalized Data 



Degrees Degrees 

Shoulder Flexion: velocity = 120 deg/sec 


Normalized Data: 1 4 Subjects Polynomial Fit to Normalized Data 



Degrees Degrees 


Figure 20. Polynomial regression fit to normalized data. 
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4.0 RESULTS AND DISCUSSION 

Fourteen subjects were measured for isolated joint torques in each axis, direction, and 
for four different velocities. This study shows how these time-consuming measure- 
ments of dynamic isolated joint torques can be predicted. A strong correlation 
between torque produced by an isolated joint and the individual's lean body mass was 
found. This relationship allows prediction of the torque-position curves for each 
isolated joint by a simple index, the lean body mass. 

To verify our regression model for predicting isolated joint torque-position curves, we 
plotted predicted versus measured torque-position curves. Figure 21 is an example of 
a prediction of an isolated joint torque-position curve (from lean body mass information 
only) plotted with measured data, the general shape of the curves of the data are 
reproduced. 


Velocity = 120 deg/sec 


Shoulder Flexion 



Figure 21. Isolated joint curve 



Degrees 

predicted versus measured. 


To get an estimate of the deviation from the measured data, a comparison computation 
was done. For each of the torque-position data sets collected, a corresponding array 
of torque values was calculated from the predicted polynomial coefficients for that 
particular data set. A torque difference vector (a difference of the actual measured 
torque value at that angle and the computed value at that point from the predicted 
polynomial coefficients) was then created, A percentage absolute value (relative to 
the maximum) of the difference array was calculated and plotted. This same analysis 
was done on the regression coefficients computed from that particular data set. The 
regression coefficients represent the ideal curve through the collected data values. 

Figures 22, 23, and 24 represent a comparison of the error of the measured data fit 
regression coefficients versus the error from the lean body mass predicted coefficients 
for one representative subject. The error of the regression coefficients reflects the 
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deviation in the fit to the measured data. In these figures, the average error of the fit to 
the measured data is 6%. For prediction of second order regression coefficients from 
lean body mass information, this is the best that one can hope to achieve. As 
expected, the predicted coefficients from lean body mass have a greater deviation 
than the regression coefficients calculated from the actual data. For the shoulder 
flexion/extension, elbow flexion/extension, and wrist flexion/extension the deviations 
are 9.8% for this individual. Figure 25 represents the same comparison for all joints 
for all individuals. When compared across the entire subject pool, the shoulder and 
elbow joint prediction equation represents about a 12% absolute value deviation. 

The wrist joint has a greater deviation. The wrist data have, on the average, one-tenth 
the range of the shoulder and elbow joints. As a result of this small range, a 
comparatively small deviation translates into a large percentage difference. The wrist 
data need further consideration. 

If only lean body mass data on a particular individual are available, our results indicate 
that the prediction equations can yield results on the order of 12% deviation from the 
actual data for all axes and directions of the elbow and shoulder joints. The results 
here could also be used to fill in missing or partial strength information for an 
individual. For instance, if only one axis of a particular joint was measured, the other 
strength information of the other axes and directions could be extrapolated from our 
equations taking into account the known information. Other correlates to dynamic 
strength prediction may exist. An area of future research is finding multiple correlates 
of strength. 



60 120 180 240 60 120 180 240 

Extension 1 Flexion 

Velocity (deg/sec) 


Figure 22. One subject's shoulder joint absolute value error of the predicted curve as 
compared to the error of the regression curve generated through the collected data. 
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Elbow Flexion /Extension 




Extension 1 Flexion 

Velocity (deg /sec) 


Figure 23. One subject's elbow joint absolute value error of the predicted curve as 
compared to the error (variation) of the regression curve generated through the 

collected data. 
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Wrist Flexion/Extension 


I Predicted from lean body mass 


■ Regression of measured data 



60 120 180 240 60 120 180 240 

Extension 1 Flexion 

Velocity (deg /sec) 


Figure 24. One subject's wrist joint absolute value error of the predicted curve as 
compared to the error (variation) of the regression curve generated through the 

collected data. 
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Figure 25. Entire subject pool. This is the absolute value error of the predicted curve 
as compared to the error (variation) of the regression curve for all subjects for all axes 

of the joints. 


Our model was based on measurements of 14 subjects. These subjects were all 
healthy young adults. Extending our lean body mass-torque regression model beyond 
this scope would require additional validation. 

Our results include several useful tables: 

1 . Tables of second order polynomial coefficients were generated from the 
measurements for each joint, axis, direction, and velocity. These tables 
provide a convenient and compact storage and retrieval mechanism to 
access our entire measured subject pool data of dynamic isolated joint 
torques (Appendix A). 

2. Tables of first order polynomial coefficients relating lean body mass of an 
individual to the average torque for an isolated joint motion were generated. 
They allow the calculation of an average torque for a particular joint, 
direction, and velocity from an individual's lean body mass (Appendix B). 
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3. Also available are second order coefficient tables characterizing the torque- 
position curves that have been normalized across our subject pool. 

These curves represent the general (normalized) trends in a torque- 
position curve and may be useful in studying torque as a function of joint 
angle in a population (Appendix C). 

To make efficient and practical use of all the coefficients, variables, and relationships 
found in this report, we have encapsulated the results of this study in a tool kit of 
software routines executing on a variety of platforms such as UNIX, DOS, and 
Macintosh machines (Appendix D). This code contains our isolated joint torque model 
(with tables of torque coefficients) that allows the prediction of dynamic isolated joint 
torques from an individual's lean body mass. 

Although our report has focused on isolated joint prediction, we did measure and 
correlate a multi-joint task (the ratchet wrench push-pull) to lean body mass. Figure 26 
indicates that multi-joint tasks may also be related to lean body mass. Percentage 
body fat is a good predictor of torque. There is a strong correlation (r > 0.92) between 
torque production capability and the lean body mass. Once a representative sample 
of a population has been measured for a particular composite motion, joint strength 
and prediction of torque capability of a particular individual may be extrapolated by 
only two measures: percentage body fat and weight. Research continues in this area. 

In addition, this ratchet data as well as all the isolated joint data are the basis of a 
graphically based human strength model being developed at the NASA-JSC Man- 
Systems Division. This model calculates end effector strength based on any arbitrary 
motions using this dynamic isolated joint information and a vector sum algorithm. 
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5.0 FUTURE DIRECTIONS 

1 . Extend isolated joint measurements to all joints. 

2. Extend the subject pool to include a more diverse group. 

3. Establish a database of dynamic isolated joint torques for general use. 

4. Examine velocity dependence in more detail. 
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7.0 APPENDICES 


Note: All data and code presented in this paper are available in 
electronic form. Contact us. 

APPENDIX A - Polynomial coefficients of angle versus torque collected 
subject. 


right_elbow 

y 

extension 


4 

60.000000 

120.000000 

180.000000 

240.000000 
flexion 

4 

60.000000 

120.000000 

180.000000 

240.000000 
right_shoulder 
x 

abduction 

4 

60.000000 

120.000000 

180.000000 

240.000000 
adduction 

4 

60.0000000 

120.000000 

180.000000 

240.000000 
flexion 

4 

60.000000 

120.000000 

180.000000 

240.000000 


2.1 1 2422E+01 
2.1 941 90E+01 
2.431055E+01 
2.676228E+01 


3.231 453E+01 
2.941 336E+01 
2.136176E+01 
3.068962E+01 


6.1 07505E+01 
4.650497E+01 
6.524007E+01 
5.786240E+01 


8.387678E+01 

8.648427E+01 

9.082626E+01 

7.684333E+01 


8.661 180E+01 
7.382412E+01 
6.001 299E+01 
5.837907E+01 


z 

lateral 

4 

60.000000 

120.000000 

180.000000 

240.000000 
medial 

4 

60.000000 
120.000000 
180.000000 
240.000000 
right_wrist 

x 


2.206026E+01 

2.191575E+01 

1.785649E+01 

1.839262E+01 


3.62771 4E+01 
4.231037E+01 
2.960975E+01 
4.6581 50E+01 


2.853215E-01 

2.263627E-01 

1.466193E-01 

-2.846204E-02 


9.382667E-02 
8.472742E-02 
2.1 1 5664E-01 
-1.243024E-01 


-2.536474E-01 

-1.961810E-02 

-5.633558E-01 

-2.265438E-01 


-1.146177E+00 
-8.942772E-01 
-9.451 21 9E-01 
-8.672966E-01 


-6.727538E-01 
-5.487776E-01 
-3.951 343E-01 
-4.092105E-01 


-1.139895E-01 
-5.408683E-02 
-5.820553E-02 
-5.881 445E-02 


8.1 51 554E-02 
4.297858E-02 
1.353765E-01 
-9.520503E-02 


-1 .952931 E-03 
-1.558132E-03 
-1.174157E-03 
7.66061 IE-05 


-1 .072887E-03 
-1 .000257E-03 
-1.517510E-03 
1 .887704E-04 


-1 .324926E-03 
-2.800043E-03 
1.123074E-03 
-2.315650E-03 


3.783874E-03 
3.240451 E-03 
3.533031 E-03 
3.563828E-03 


1.768133E-03 
1 .41 9674E-03 
1.065975E-03 
1. 30531 6E-03 


-1.104627E-03 
-1 .786900E-03 
-7.560909E-04 
5.387657E-04 


-4.386652E-04 
-1.839004E-03 
-1 .624337E-03 
-1. 27341 2E-03 


0.000000E+00 

0.000000E+00 

O.OOOOOOE+OO 

0.000000E+00 


0.000000E+00 

O.OOOOOOE+OO 

0.000000E+00 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


data for one 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 
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radial 

4 

60.000000 

120.000000 

180.000000 

240.000000 
ulnar 

4 

60.000000 
120.000000 
180.000000 

240.000000 

y 

extension 

4 

60.000000 
120.000000 
180.000000 

240.000000 
flexion 

4 

60.000000 
120.000000 
180.000000 

240.000000 
z 

pronation 

4 

60.000000 
120.000000 
180.000000 

240.000000 
supination 

4 

60.000000 
120.000000 
180.000000 
240.000000 


1.141442E+01 

1.020690E+01 

9.970860E+00 

8.380389E+00 


1.137878E+01 

1.174232E+01 

1.070129E+01 

8.977362E+00 


5.7371 69E+00 
5.026070E+00 
4.333305E+00 
4.2201 27E+00 


1. 27601 3E+01 
1.155679E+01 
1.146068E+01 
1.0581 81 E+01 


7.31 1457E+00 
7.3361 91 E+00 
6.619156E+00 
6.05881 6E+00 


5.309766E+00 

5.978102E+00 

5.724769E+00 

6.056996E+00 


9.410274E-04 
-5.753472E-03 
-6.447901 E-02 
-7.270560E-02 


9.434845E-02 

1.582629E-01 

1.767397E-01 

1.348160E-01 


7.455280E-02 
5.843269E-02 
7.8331 46E-02 
8.03501 9E-02 


2.093285E-03 
-2.462766E-02 
-1.606610E-02 
-5.0641 55E-02 


-3.482732E-02 
-4.6641 70E-02 
-2.951 976E-02 
-2.663488E-02 


4.141473E-02 
3.62561 9E-02 
3.90581 2E-02 
4.260557E-02 


-3.824661 E-03 
-3.008693E-03 
-1.006891 E-03 
-3.003046E-04 


-2.236507E-03 
-3.33841 5E-03 
-2.987377E-03 
-1 .860539E-03 


-3.73401 9E-04 
-2.478556E-04 
3.712687E-04 
-1. 00261 0E-04 


-1.65355 IE-03 
-2.024493E-03 
-1.736379E-03 
-1.466442E-03 


-6.820944E-04 

-6.950588E-04 

3.903088E-04 

-3.767824E-04 


-3.062696E-04 

-5.193476E-04 

-3.756247E-04 

-4.954465E-04 


0.000000E+00 

O.OOOOOOE+OO 

0.000000E+00 

0.000000E+00 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 
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APPENDIX B - Coefficient table relating lean body mass to average torque for a particular velocity. 


right_elbow 

y 

extension 
4 

60.0000 
120.0000 
180.0000 

240.0000 
flexion 
4 

60.0000 

120.0000 
180.0000 

240.0000 
right_shoulder 
x 

abduction 
4 

60.0000 

120.0000 
180.0000 

240.0000 
adduction 
4 

60.0000 

120.0000 
180.0000 

240.0000 

y 

extension 

4 

60.0000 

120.0000 
180.0000 

240.0000 
flexion 
4 

60.0000 

120.0000 
180.0000 

240.0000 
z 

lateral 
4 

60.0000 

120.0000 
180.0000 

240.0000 
medial 
4 

60.0000 

120.0000 
180.0000 
240.0000 
right_wrist 
x 


-1 .78031 1E+01 
-1.777659E+01 
-1 .575342E+01 
-1.138565E+01 


-1 .1701 13E+01 
-1.499996E+01 
-1.190419E+01 
-9.734751 E+00 


-1.162713E+01 
-1 .00581 9E+01 
-9.859225E+00 
-8.557056E+00 


-1 .840421 E+01 
-1.617023E+01 
-2.081 728E+01 
-1.698798E+01 


-2.438441 E+01 
-2.267868E+01 
-1.129235E+01 
-2.071 749E+01 


-1.827386E+01 
-2.1 1 8529E+01 
-9.36631 6E+00 
-1.645339E+01 


-8.93631 5E+00 
-7.1 74304E+00 
-4.212536E+00 
-6.357865E+00 


-1.540060E+01 
-1 .11 4382E+01 
-1.053552E+01 
-1 .421 859E+01 


6.688885E-01 
6.585007E-01 
6.0431 51 E-01 
5.154829E-01 


6.21 4280E-01 
6.524455E-01 
5.64501 5E-01 
4.945650E-01 


6.52501 9E-01 
5.898796E-01 
5.921 002E-01 
5.585406E-01 


8.565667E-01 
8.1 01 834E-01 
8.792474E-01 
8.091 762E-01 


9.66741 0E-01 
9.26651 8E-01 
6.457043E-01 
8.310267E-01 


8.346249E-01 
8.891 425E-01 
6.267306E-01 
7.376269E-01 


4.264861 E-01 
3.740323E-01 
3.108540E-01 
3.512995E-01 


6.548921 E-01 
5.727975E-01 
5.327551 E-01 
6.126382E-01 


0.000000E+00 

0.000000E+00 

0.000000E+00 

O.OOOOOOE+OO 


0.000000E+00 

0.000000E+00 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.GOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 
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radial 

4 

60.0000 

120.0000 

180.0000 

240.0000 
ulnar 

4 

60.0000 

120.0000 
180.0000 

240.0000 

y 

extension 

4 

60.0000 

120.0000 
180.0000 

240.0000 
flexion 

4 

60.0000 

120.0000 
180.0000 

240.0000 
z 

pronation 

4 

60.0000 

120.0000 
180.0000 

240.0000 
supination 
4 

60.0000 

120.0000 
180.0000 
240.0000 


-1.268475E+00 
-1.005637E+00 
-1.270290E+00 
-1. 77091 2E+00 


-2.900480E+00 

-3.978939E+00 

-2.994360E+00 

-3.975077E+00 


-1.127101 E+00 
-7.367646E-01 
-1 .423397E+00 
-1.275451 E+00 


4.74761 6E+00 
-4.260902E+00 
-2.04081 7E+00 
-3.6541 24E+00 


-6.592008E-01 

-1.366900E+00 

-1.853492E+00 

-2.182193E+00 


-2.369235E+00 
-2.446581 E+00 
-2.280463E+00 
-2.641 539E+00 


1.357690E-01 
1.208873E-01 
1.214390E-01 
1 .26391 8E-01 


1 .9329001-01 
2.042846E-01 
1.823847E-01 
1.952298E-01 


8.416442E-02 
7.451 71 5E-02 
8.59421 3E-02 
8.026770E-02 


2.070452E-01 
1 .903496E-01 
1.303384E-01 
1.715256E-01 


7.739327E-02 

9.181561E-02 

1.013228E-01 

1.070073E-01 


1 .065645E-01 
1 .08827 9E-01 
1 . 070542 E-01 
1 .1 19357E-01 


0.000000E+00 

0.000000E+00 

0.000000E+00 

O.OOOOOOE+OO 


0.000000E+00 

0.000000E+0O 

O.OOOOOOE+OO 

0.000000E+00 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO i 

O.OOOOOOE+OO = 

O.OOOOOOE+OO 1 

O.OOOOOOE+OO 


0 OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 
0^000000 E+00 
O.OOOOOOE+OO 
O.OOOOOOE+OO 


40 


'I III II !| IINIilPlinil I I I Ill II II III 1 III ill! 111 II Ml !l IHI'IIIIIMII » II 



Normalized torque coefficients for entire subject pool. 


Appendix C - 
right_elbow 

y 

extension 

4 

60.0000 

120.0000 

180.0000 

240.0000 
flexion 

4 

60.0000 

1 20.0000 

180.0000 

240.0000 
right_shoulder 
x 

abduction 

4 

60.0000 

120.0000 
180.0000 

240.0000 
adduction 
4 

60.0000 

120.0000 
180.0000 

240.0000 

y 

extension 

4 

60.0000 

120.0000 
1 80.0000 

240.0000 
flexion 

4 

60.0000 

120.0000 
1 80.0000 

240.0000 
z 

lateral 

4 

60.0000 

120.0000 
180.0000 

240.0000 
medial 

4 

60.0000 

120.0000 
180.0000 
240.0000 
right_wrist 
x 

radial 


2.382654E-01 

3.998269E-01 

6.186674E-01 

6.797335E-01 


9.436899E-01 
9.559348E-01 
1.024970E+00 
1.1811 90E+00 


1 .63821 7E+00 
1.548579E+00 
1.754382E+00 
1.667307E+00 


8.051 392E-01 
8.1 1 4389E-01 
8.734385E-01 
7.575982E-01 


6.296859E-01 

4.782509E-01 

7.084687E-01 

6.559435E-01 


1.873006E+00 
1.857643E+00 
1.801 71 9E+00 
1 .71 0937E+00 


1 .061 949E+00 
1 .041 042E+00 
1.007360E+00 
9.669892E-01 


1.022044E+00 
1.041949E+00 
1.029437E+00 
1.021 175E+00 


1 .609246E-02 
1.352749E-02 
7.465366E-03 
5.123406E-03 


1 .385255E-03 
2.457487E-03 
7.60521 7E-04 
-2.677844E-03 


-9.23421 3E-03 
-6.3281 32E-03 
-1.159083E-02 
-9.909423E-03 


4.159272E-03 
4.182187E-03 
1 .71 7887E-03 
4.319388E-03 


7.588339E-03 
9.29431 4E-03 
4.127454E-03 
3. 1 29994 E -03 


-1 .086977E-02 
-1.102275E-02 
-1.150634E-02 
-9.788573E-03 


-5.109649E-03 
-4.551 276E-03 
-4.989951 E-03 
-4.986339E-03 


4.982772E-03 

4.882844E-03 

4.798532E-03 

3.547545E-03 


-7.057701 E-05 
-6.401 985E-05 
-2.948732E-05 
-1 .389026E-05 


-7.300826E-06 
-1.990186E-05 
-1.120966E-05 
3.1 23384E-06 


1 .584508E-05 
-3.656676E-06 
2.1 8831 2E-05 
1 .595937E-05 


-1.923175E-05 

-2.022592E-05 

-2.843898E-06 

-1.588424E-05 


-3.253375E-05 
-3.527058E-05 
-1.1471 16E-05 
-1.1 57555 E-06 


2.096022E-05 

2.299863E-05 

2.964644E-05 

2.236658E-05 


-1 .644733E-05 
-1 .793456E-05 
9.01 7851 E-06 
2.867537E-05 


-3.170801 E-05 
-5.522752E-05 
-4.722543E-05 
-3.61 5252E-05 


0.000000E+00 

0.000000E+00 

0.000000E+00 

O.OOOOOOE+OO 


0.000000E+00 

O.OOOOOOE+OO 

0.000000E+00 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 
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4 

60.0000 

120.0000 

180.0000 

240.0000 
ulnar 

4 

60.0000 

120.0000 
180.0000 

240.0000 

y 

extension 

4 

60.0000 

120.0000 
180.0000 

240.0000 
flexion 

4 

60.0000 

120.0000 
180.0000 

240.0000 
z 

pronation 

4 

60.0000 

120.0000 
180.0000 

240.0000 
supination 
4 

60.0000 

120.0000 
180.0000 
240.0000 


1.181802E+00 

1.107246E+00 

1.036549E+00 

1.039976E+00 


-1. 67591 3E-02 
-1 .877522E-02 
-1.984581 E-02 
-2.084383E-02 


1.066880E+00 
1.087365E+00 
1.1 14945E+00 
1.07331 1E+00 


9.806204E-01 
9.499655E-01 
9.1 77343E-01 
9.216082E-01 


1.084562E+00 
1.109554E+00 
1.091 251 E+00 
1.088930E+00 


1.01 0391 E+00 
9.950561 E-01 
1.009005E+00 
9.778854E-01 


1.024255E+00 
9.929054E-01 
1 .009055E+00 
9.361 496E-01 


9.789578E-03 
1 .1 14404E-02 
1.416128E-02 
1.608139E-02 


8.9601 70E-03 
9.603754E-03 
8.428668E-03 
9.430406E-03 


4.1 2461 7E-04 
5.107403E-04 
-5.016060E-05 
-9.37731 2E-04 


5.380042E-03 
5.7331 34E-03 
5.072503E-03 
6.124177E-03 


-7.277947E-03 
-7.556385E-03 
-7.271 954E-03 
-6.8241 19E-03 


-6.016592E-04 

-4.386299E-04 

-3.919646E-04 

-3.424259E-04 


-3.1 74051 E-04 
-3.558657E-04 
-3.47691 9E-04 
-3.90451 5E-04 


-5.493223E-05 
-4.643842E-05 
-3.641 872E-06 
-2.801 543E-05 


-1.078120E-04 
-1 .394872E-04 
-1 .093535E-04 
-1 .0801 13E-04 


3.692308E-06 
7.248831 E-06 
-6.918676E-06 
1.015605E-05 


-3.147846E-05 
-3.043222E-05 
-5.40331 2E-05 
-1.593118E-05 


O.OOOOOOE+OO 

0.000000E+00 

0.000000E+00 

0.000000E+00 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 
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//*add to this list as more joints become available and increment JC the joint count*/ 
static char joint_name[MAXJOINT|[80] = {"right_shoulder", 

"nghtjslbow", 

"right_wrist" }; 

static AxisTable AxisMapfl = 

{ 

"ABDVEL", "abduction", "x", 1, 

"ADD VEL", "adduction", V,-1, 

"FLEXVEL", "flexion", "y", 1, 

"EXT VEL", "extension", "y*,-1, 

"EXT ROT, "lateral", "z", 1, 

"INTROT, "medial", "z", -1, 

"PRONATION VELOCITY", "pronation", "z", 1, 

"SUPINATION VELOCITY", "supination", "z", -1, 

"RADIAL DEV", "radial", V, 1, 

•ULNAR DEV", "ulnar", V,-1, 

0 , 0 , 0 , 0 

}; 

int AxisTableSize = 10; 
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CalculatelndividualCoeffTable (itable, ptable, Itable, bf, weight, count) 


CoeffTable itableQ; 
CoeffTable ptableQ; 
CoeffTable ItableQ; 
float bf; 
float weight; 
int count; 


Function: 

Given the population normalized table (ptable) and the 
lean body mass to torque relationship table ( Itable ), this program will 
calculate the the itable ( individual person's table) from his/her 
body fat and weight. 

The itable is calculated by computing the average torque for each 
joint direction and velocity from the Itable. The itable (lean body mass 
table) contains coefficients which describe the relationship between 
the average torque and the lean body mass. To use this table, first 
a lean body mass is computed from the weight and % body fat This value is 
plugged into the coefficients looked up in the Itable. This yields the 
average torque for that joint direction and velocity. After this torque is 
computed, the normalized torque coefficients for that joint, direction and 
velocity are looked up in the ptable and multiplied by the average torque 
computed from the Itable. These coefficients are then stored into the itable 
and returned. 

Parameters: 
itable - individuals table 
ptable - population/ normalized table 
Itable - lean body mass to average torque table 
bf - body fat. 
mass - mass of the person. 


{ 

int axis; 
ait dir; 

int eq_count; 
int i, j , z; 

float lean_body_mass, average; 
float ComputeCoeff 0 

/*lean body mass calculation*/ 
lean_body_mass = weight * ( 1 .0 - bf );; 
printf ( "Coefficients calculated for lean body mass %f , lean_body_mass); 
/*for each joint*/ 
for (j = 0 ; j < count; j++) 

{ : .. 

/* copy the name and axis count to the itable*/ 
strcpy ( itable[fl.joint_name, ptable[j].joint_name); 
itableQ], axis_count = ptabieQ].axis_count; 

/* for each axis in this joint*/ 

for( axis = 0; axis < ptableQ]. axis_count; axis++ ) 

{ 

/*copy the axis number*/ 
itableQ].axis[axis] = ptableQ].axis[axis]; 
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/"for each direction in this axis*/ 
for( dir = 0; dir < 2; dir++) 

/* copy the direction name and equation count to the ptable*/ 
strcpy ( itable[j].function[axls][dir].direction_name p 
ptable[j].1ijnction[axis][dir].direction_name); 
itable[j].function[axis][dir].eq_count = 
ptable[j].function[axis][dir].eq_count; 
eq_count = ptableO]function[axis][dir].eq_count; 

/*for each regression equation compute the individuals eqns*/ 

tbr( i = 0; i < eq_count; I++ ) 

/*copy the velocity information*/ 
itable[j].funct'on[axis][dir].velocity[i] = 
ptable[j].function[axis][dir].velocityO]; 

/*compute the average for torque for this dir, velocity*/ 

/*ffom the lean body mass information and the Itable */ 
average = 

ComputeCoeff (ltable[j].function[axis][dir].coeff[i], 
leanjx>dy_rnass); 

/•multiply the coefficients of the pop table by average*/ 

/*and store in the itable */ 

MultiplyCoeff (itableffl .fu notion [axi s] [dir], coeffp} , 
ptabie[0.function[axis][dir].coeff[i], 
average); 

} Tend fori*/ 

}/*endfbrdir*/ 

y* end for axis*/ 

}/*end for j-joint*/ 


} 

WriteCoeffTableToFile( table, ext, count ) 
CoeffTable table []; 
char extfl; 
int count; 


Function: 

Given a table pointer a file extension, and the number of joints, it will write 
the coefficients of the table specified into files with the names 
of the files taken form the joint names in joint_name table with the extension 
in variable ext. 

This is an example of the file format: 

right_elbow 

Y 

extension 

4 

60.000000 1.743286e+01 1.177417e+00 -5.163818e-03 0.000000e+00 

120.000000 2.869938e+01 9.709965e-01 -4.5953136-03 0.000000e+00 

180.000000 4.1 10020e+01 4.959499e-01 -1 ,958944e-03 0.000000e+00 

240.000000 3.991 397e+01 3.008465e-01 -8.156364e-04 0.000000e+00 
extension 
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4 

60.000000 1.743286e+01 1.177417e+00 -5.163818e-03 0.000000e+00 

120.000000 2.869938e+01 9.709965e-01 -4.59531 3e-03 0.000000e+00 

180.000000 4.1100206+01 4.959499e-01 -1.958944e-03 0.000000e+00 

240.000000 3.9913976+01 3.0084656-01 -8.1563646-04 0.0000006+00 

line 1 : joint name 
line 2: axis 
line 3: direction 

line 4: number of regression equations, 
line 5: velocity, coefficient 1 coefficient 2 ... 
this repeats until the entire joint is described. 

Parameters: 

table - storage table for the coefficient values. 

ext - file extension the data to be stored in ( eg ffc , ntc or Ibc ) 

count - number of joints. 


FILE *fp; 

char filename[80]; 
char line[132]; 
int I, j ; 
int axis; 

char axisname[32]; 
char dimame[80j; 
int ano; 
int dir; 

Int eq_count; 
int joint; 

/* for each joint in the joint names, 

write the file for that joint 

for the table specified.*/ 

for ( joint = 0 ; joint < count ; joint++) 

/*open a file for this joint using names form the 
joint_name table 
7 


strcpy ( filename, joint_name[joint]); 
sprlntf( filename, "%s.%s“, filename, ext ); 
fp = fopen( filename, “w" ); 
if(fp == NULL) 

{ 

printf( "failed open on %s\n", filename ); 
return ( 0 ); 

} 

/* write out name of joint */ 
fprintf( fp, "%s\n\ tableQoint].joint_name ); 
for (axis = 0 ; axis < table[joint].axis_count ; axis++ ) 


if ( table[joint].axis[axis]== 0 ) fprintf(fp, "X\n“); 
if ( tablepnt].axis[axis]== 1 ) fprintf(fp, "V\n"); 
if ( table[jointj.axis[axis]== 2 ) 1printf(fp, "Z\n"); 


/* direction name */ 
1br( dir = 0; dir < 2; dir++ ) 

{ 
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fprintf(fp, "%s\n", tableDoint].function[axis][dir].direction_name ); 
I* number of force equations for this direction */ 


fprintf( fp, "%d\n", 

table[joint] .function[ax!s] [dir].eq_count ) ; 
eq_count = table[joint].function[axis][dir].eq_count; 


} 


} 


} 

} /* axis loop */ 

fdose(fp); 

} /*end for each joint*/ 


for( i = 0; I < eq_count; i++ ) 

fprintf( fp, "%f %e %e %e %e \n", 

tableOointj.function [axis] [d ir] .velocity [i] , 
tableQoint] .function [axis] [dir] .coefffi] [0] , 
table[joint].function[axisi[dir].coeff[i][1]. 
table[joint].function[axis][dir].coeff[i][2], 
tableyoint] .function [axis] [dir] .coefffi] [3] ) ; 


TorqueSummary ( table, bf, weight, count) 

CoeffTable tableO; 

float bf; 

float weight; 

int count; 



Function: Print out a quick-look table of the average torque values 
for each of the joints in the joint chain. This function takes the lean 
body mass table coefficients the %bf, and weight of an individual and computes 
the average torque for joint, direction. All velocities are averaged before 
being printed. This means that the dynamic information is not printed in order 
to keep it a quick look table. 

Parameters: 

table - the lean body mass table, 
bf -%bodyfat(0-1) 
weight - weight of the person, 
count - joint count. 


{ 

int axis; 
int dir; 

int eq_count; 
int i.j; 

float average, avg; 
float ComputeCoeff 0; 
float lean_body_mass; 
lean_body_mass = weight * ( 1 .0 - bf ); 
printf (“\n\n"); 

printf (■ W); 

printf (" Weight %f Body Fat %f Lean Mass %1\n', weight, bf, lean_body_mass); 
for (j = 0 ; j < count; j++) 

printf( "\n\nJoint Name: %s\n", table[j].joint_name ); 
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for( axis = 0; axis < table[j].axis_count; axls++ ) 
for( dir = 0; dir < 2; dir++ ) 

{ 

printf( "Direction: %s", 

table[j].function[axis][dir].direction_name); 

eq_count = tablep].function[axis][dir].eq_count; 
for( i = 0; i < eq_count; i++ ) 

{ 

/•compute the average for torque for this dir, velocity*/ 
average = 

ComputeCoeff (table|j].function[axis][dir].coeff[i], 
lean_body_mass); 
avg +=average; 


} 

} 


} 


} 

printf ( " Average Torque: %f\n", avg / eq_count); 
avg = 0.0; 


printf ("■ 

} 


An"); 


ReadCoeffTable( table, ext, count ) 


Function: 

Given a table pointer a file extension, and the number of joints, it will fill 
the coefficient table specified with the values in the file according to 
the joint names in the joint_name table defined. 


Parameters: 

table - storage table for the coefficient values, 
ext - file extension the data are stored in (eg ffc , ntc or Ibc ) 
count - number of joints. 


CoeffTable tableQ; 
char extQ; 
int count; 

{ 

FILE *fp; 

char filename[80j; 
char line[132j; 
int i, j ; 
char *ier; 
int axis; 

char axisname[32j; 
char dirname[80j; 
int ano; 


int dir; 

int eq_count; 


int joint; 


/* initialize force table 7 
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lnitJointCoeffTable( table ); 

/* for each joint in the joint names, 
load the file for that joint 
into the table specified.*/ 

for ( joint = 0 ; joint < count ; joint++) 

' strcpy (filename, joint_nameDointj); 

sprintf( filename, *%s.%s*, filename, ext ); 
fp = fopen( filename, V ); 
rf( fp == NULL) 

^ printf( "failed open on %s\n", filename ); 
return(O); 


r read in name of joint 7 

ier = fgets( line, 132, fp ); 

if( ier == NULL ) return( 0 ); 

sscanf( line, "%s\n", tableDoint].joint_name ); 

/* read remaining data in file 7 
axis = 0; 
while( 1 ) 

ier = fgets( line, 132, fp ); 
if( ier == NULL ) break; 

/* axis number 7 

sscanf( line, "%s", axisname ); 

ano = FTGetAxisNumber( axisname[0] ); 

if( ano == -1 ) 

^ printf( "invalid axis name %s\n", axisname ); 
printf( "using default X axis\n" ); 
ano = 0; 


table[joint].axis[axis] = ano; 
fbr(j = 0;j<2; j++) 

ier = fgets( line, 132, fp ); 
if( ier == NULL ) break; 
f* direction name 7 
sscanf( line, '%s", dirname ); 

/* get direction code 0 = negative, 1 = positive, 2 = error / 
dir = GetAxisDirection( dimame ); 
if( dir == 2 ) dir = 0; 

strcpy( table[joint].function[axis][dir].direction_name, 
dimame); 

ier = fgets( line, 132, fp ); 
jf( ier == NULL ) break; 


/* number of force equations for this direction 7 
sscanf( line, "%d", 

&tableOoint].function[axis][dir].eq_count ); 
eq_count = table[joint].function[axis][dir].eq_count; 
for( i = 0; i < eq_count; I++ ) 

ier = fgets( line, 132, fp ); 

if( ier == NULL ) break; /* break out of for loop 7 
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sscanf( line, "%f %e %e %e %e", 

&tableQoint] .function [axis] [dir] .velocity [i] , 
&tablepoint].function[axis][dir].coeff[i][0], 
&tablepointj .function [axis] [dir] .coefffi j[1 j , 
&table0oint].function[axis][dirj.coeff[i][2], 
&tablepoint].function[axis][dir].coeff[i][3] ); 


} 

/* break out of the while loop */ 
if( ier == NULL ) break; 
axls++; 

tablepointj.axis count = axis; 

> 

fclose (fp); 

} /*end for each joint*/ 


} 

int InitJointCoeffTable ( table ) 

Function: Given a pointer to the coefficient table, this function 
will put all the appropriate initial values into the table. 
Parameters: 

table - pointer to the coefficient table to be initialized. 

CoeffTable tabled; 

{ 

int i, j, k, I, m ; 

/*for each joint int the table*/ 
for (m = 0; m < MAXJOINT; m++) 

{ 

table[m].jolnt_name[0] = 0; 
table[m].axls_count * 0; 
table[m].axls[0] = 0; 

table[mj.axis[lj =1; 

table[m].axls[2] = 2; 

/*for each axis in that joint*/ 
fbr( i = 0; I < AXIS; I++) 

{ 


/*for each direction in that axis*/ 

■fbr( j = 0; j < DIR; j++ ) 

{ - 

table[m].function[i]p'].eq_count = 0; 
table[mj.function[i]p].direction_name[0] = 0; 
/* for each velocity in that direction*/ 
for( k = 0; k < VELOCITY; k++ ) 

{ 


table[m] .function [ij pj .velocity [k] = 0.0; 
for( I = 0; I < COEFF; I++ ) 

table[m], function p]p].coeff[k][l] = 0.0; 

} 


} 


} 


} 


} 
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UstCoeffTable( table, label, count ) 
CoeffTable tabled; 
char labelQ; 
int count; 


Function: Given a pointer to a table, this function 
writes the function to stdout. 

Parameters: 

table - pointer to a table. 


* * **** * *********** * A A AAAA ****** ** 


{ 

int axis; 

int dir; 

int eq_count; 

int i,j; 

printf ( "%s\n", label); 

for (j = 0 ; j < count; j++) 

{ 

printf( "Joint Name: %s\n", table[j].joint_name ); 

printf( "DOFs : %d\n", tablefi].axis_count ); 

for( axis = 0; axis < tableQ].axis_count; axls++ ) 

{ 

fbr( dir = 0; dir < 2; dir++ ) 

{ 

printf(" Direction: %s\n", 

table[j] .function [axis] [dir] ,direction_name ) ; 
printf( " Number of Functions: %d\n", 
table[j].function[axis][dir].eq_count ); 

printf(" Velodty Coefficients — -\n"); 

eq_count = table[j].function[axis][dir].eq_count; 
for( i = 0; i < eq_count; i++ ) 

{ 

printf(" %f %E %E %E %E\n", 
table[]].function[axis][dir].velocity[i], 
table|j] .function [axis] [dir] .coeff[i] [0] , 
tableQ].function[axis][dir].coeff[i][1 j, 
tabled .function [axis] [dirj .coeffp] [2] , 
table0.function[axis][dir].coeffp][3] ); 

} 

} 

> 

} 

} 


float ipower( x, y ) 

float x; 

inty; 


**************** 


float ipower( x, y ) - compute x to the y power 
Function: 

Compute the value of x raised to the yth power. 


{ 

float z; 
z = 1 .0; 

while(y-){z*=x;} 

retum(z); 

} 
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MultiplyCoeff ( coeffl , coeff2, value) 
float coeffl D; 
float coeff2Q; 
float value; 

Function: 

Multiply the array of coeffs by the value specified, 
coeffl = coeff2 * value; 

Parameters: 

coeff - an array of float values that are the coefficients to 
an n degree polynomial, 
value - the value to multiply with. 


{ 

inti; 

for( i = 0; i < COEFF; i++ ) 

{ 

coeffl [i] =coeff2[i] * value; 

} 

} 

float ComputeCoeff ( coeff, value ) 
float coeff Q; 
float value; 


Function: 

Solve the coefficients for the value specified. 
Parameters: 

coeff - an array of float values that are the coefficients to 
an n degree polynomial. 

val ue - th e x value of the function to compute the y. 

{ 

inti; 

float answer; 

answer = coeff [0]; 

fbr( i = 1 ; i < COEFF; i++ ) 

{ 

answer += ipower( value , I ) * coeffp]; 

} 

return ( answer); 


Int FTGetAxisNumber( name ) 
char name; 

int FTGetAxisNumber( name ) Function: 

Given an axis name return axis number associated with name. 


char a; 

a = toupper( name ); 
if ( a == 'X' ) retum( 0); 
if ( a == V ) return ( 1 ); 
if(a=='Z) retum(2); 
retum( -1 ); 
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} 


n f. 1 ft i i . lllllMmii»HH»m» < » t Ht * *»»** * » **»»» > » ** ***** 

Function: Given the name of the direction, return a direction number. 

Parameters: name - direction name. 


char *name; 

{ 

int i; 
intk; 
i = 0; 

while( AxisMap[i].input_name[0] ) 

k = strlen( AxisMap[i].output_name ); 

if( strncmp( name, AxisMap[i].output_name, k ) == 0 ) 

* if( AxisMapp] .direction < 0 ) return( AxisMapp]. direction + 1 ); 
else retum( AxisMappJ.direction ); 

} 

i++; 

} 

retum(2); 

} 



GetVelocityCode( joint, axis, dir, velocity ) 

Function: 

Get the code number associated with a given velocity 
found in the force function velocity table. The code 
number is selected by comparing the given velocity with 
the velocities in the table and chosing the index number 


of the closest value. 


int GetVelocityCode( force_table, joint, axis, dir, velocity ) 

CoeffTable force_table[]; 

int joint; 

int axis; 

int dir; 

float velocity; 

{ 

int i; 

int count; 
float mid; 

i I* 

size of the velocity is the same as the number of 
equations within the given axis and direction 
7 

count = force_table[joint].function[axis][dir].eq_count; 
p check for special case of zero vel. 7 
/*no equations7 
if (count == 0 ) return ( -1 ); 

if( velocity < force_table[joint].function[axis][dir].velocity[0] ) 
retum(O); 

/* scan table 7 
for( i = 0; i < count - 1 ; i++ ) 

* if( (velocity > force_tableUoint].function[axis][dir].velocityp]) && 
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{ 

printf ( " \n\n"); 

printf ( "Fat Menu ===========> :-) \n"); 

printf ( "0 -> exit. \n"); 

printf ( “1 -> Input a new bf and weighAn\n“); 

printf ( "2 -> Print Population Coefficients^"); 

printf ( "3 -> Print Lean Body Mass coefficients^"); 

printf ( "4 -> Print This individuals coefficients\n\n"); 

printf ( "5 -> Write out (*.ffc files) individuals coefficients^"); 

printf ( "6 -> Print out a torque summary for this individualAn"); 
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} 

main (argc, argv) 
int argc; 
char *argvO; 

/ 

Function: main program. Menu selection and execution. 

Parameters: body fat and weight. 

AAAAAAAAAAAAAAAAAAAAAAAAAAA * ****** 1 " 1 " 1 

|AAAAAAAAAA"* **- AAAA r A *** *AA»AAAA*AAA»AAAAAAA»* ** AAA AAA** * *** * * 

type definitions for CoeffTable data structure see lean.h. 

Decription of force table data structure. 
table[joint] .function [axis] [di r] 
force functions associated with a joint, these 
are indexed by the axis number and the direction 
code. 

table[joint].joint_name 
name of joint 80 characters max. 
table[joint].axis[axis] 

table of axis numbers (DOFs) for the joint. 
table[joint].axis_count 
number of axes or DOFs for the joint, 
function [axis] [di r] . coeff[vid] [coeff] 
force function coeffcients index by the velocity code (vid) 

(see velocity code) and coefficient index (1-4). 
function[axis][dir].velocity[vid] 
velocity value associated with a given velocity code (vid). 
function [axis] [dir]. eq_count 
number of force function equations associated with 
a given axis and direction, 
function [axis] [dir] . directionjiame 

name of a given direction, i.e. extension, flexion, etc. 

^ !>< ll li *< i iU*>**>*«»*****»»* »»*** *»*»**»»l>*«*«*«*«*********** ' * **** ******** A *** ********* t ************* < ‘ J ' A ^ 

CoeffTable pop_table[MAXJOINT|; /*normalized population coefficients*/ 

CoeffTable lean table[MAXJOINT]; ^relationship: mean torque to lean mass*/ 

CoeffTable ind_table[MAXJOINTj; /*Calculated coeff table based on 
weight and %body fat*/ 

float bf; /* body fat percentage*/ 

float weight; /*weight */ 

int done; /loop control */ 

int option; /"user input */ 

#if (PC || UNIX) 

/* process command line arguments*/ 
if ( argc < 3 ) 

^ printf ( "usage: lean <%body_fat> <weight> \n"); 
exit( 0 ); 

} 

#endif 
#if MAC 
printf ( "\n\n"); 

printf ("input body fat (0-1) and <weight(kg):"); 
scanf ( "%f %f , &bf, &weight); 
printf ("%f%f,bf, weight); 

#else 

bf = atof ( argv[1] ); 
weight = (float) atof ( argv[2]); 
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#endif 
if (bf > 1) 

{ 

printf ( " invalid % body fat :-) enter value between 0-1 \n"); 
exit( 0 ); 

} 

/*read the two tables 1 . Population table 2. lean body table*/ 
ReadCoeffTable ( pop_table , "ntc" , JC ); 

ReadCoeffTable ( lean_table, "Ibc", JC ); 

/“calculate this individual's table*/ 

InitJointCoeffTable ( ind_table ); 

CalculatelndividualCoeffT able 
(ind_table, pop_table, lean_table, 
bf, weight, JC );” 
done = FALSE; 
while ( Idone ) 

{ 

menu 0; 

scanf ( "%d", &option); 
switch (option) 

{ 

case 0: /*ex it*/ 

printf ("May the Vlean forceV be with you. \n"); 
done = TRUE; 
break; 

case 1 : /*lnput new values*/ 

printf ("Input body fat (0-1) and weight:"); 
scanf ("%f %f, &bf, &weight); 

Calc ulate lndividualCoeffTable 
(ind_table, pop_table, lean_table, 
bf, weight JC ); 
break; 

case 2: /*Print Population Coefficients*/ 

ListCoeffTable ( pop_table, "Population Table", JC); 

break; 

case 3: /*Print Lean Body mass coefficients */ 

ListCoeffTable ( lean_table, "Lean Body mass -> average torque", JC); 
break;; 

case 4: /*Print Individuals coefficients*/ 

ListCoeffTable ( ind_table, "Individuals coefficients", JC); 
break; 

case 5: /*Write out individuals coefficient file*/ 

WriteCoeffTableToFile( ind_table, "ffc", JC); 
printf (" *.ffc files are written\n“); 
break; 

case 6: /“Write out a summary table for this individual*/ 

TorqueSummary ( lean_table, bf, weight, JC); 
break; 
default: 

printf (" TRY AGAIN, (enter options 0 ->6) \n“); 
break; 

} 

} 

} 


56 




REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden for this collection of Information is estimated to average 1 hour per response, including the time for reviewing Instructions, searching existing data sources, gathering and 
maintaining the data needed, and completing and reviewing the collection of information Send comments regarding this burden estimate or any other aspect of this collection of information, 
including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 
22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (Q7Q4-Q IBB), Washington, DC 20503- 

1. AGENCY USE ONLY (Leave blank) 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED ' 

June 1992 Technical Paper 

4 TITLE AND SUBTITLE 

Correlation and Prediction of Dynamic Human Isolated 
Joint Strength From Lean Body Mass 

6. AUTHOR(S) 

Abhilash K. Pandya, Scott M. Hasson, Ann M. Aldridge, 

James C. Mai da, and Barbara J. Wool ford 


11. SUPPLEMENTARY NOTES 

Abhilash K. Pandya and Ann M. Aldridge; Lockheed Engineering & Sciences Company, 
Houston, Texas; Scott M. Hasson: Texas Women's University, Denton, Texas; James C. 

Maida and Barbara J. Woolford: Lyndon B. Johnson Space Center, Houston, Texas. 


13 ABSTRACT (Maximum 200 words) 

A relationship between a person's lean body mass and the amount of maximum torque that 
can be produced with each Isolated joint of the upper extremity was Investigated. The 
maximum dynamic Isolated joint torque (upper extremity) on 14 subjects was collected 
using a dynamometer multi-joint testing unit. These data were reduced to a table of 
coefficients of second degree polynomials, computed using a least squares regression 
method. All the coefficients were then organized Into look-up tables, a compact and 
convenient storage/retrieval mechanism for the data set. Data from each joint, 
direction and velocity, were normalized with respect to that joint's average and merged 
Into files (one for each curve for a particular joint). Regression was performed on 
each one of these files to derive a table of normalized population curve coefficients 
for each joint axis, direction and velocity. In addition, a regression table which 
Included all upper extremity joints was built which related average torque to lean body 
mass for an individual. These two tables are the basis of the regression model which 
allows the prediction of dynamic Isolated joint torques from an Individual's lean body 
mass. 


12b. DISTRIBUTION CODE 


12a. DISTRIBUTION / AVAILABILITY STATEMENT 

Unclassified - Unlimited 
Subject Category 54 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 
Man-Systems Division 

National Aeronautics and Space Administration 
Lyndon B. Johnson Space Center 
Houston, Texas 77058 

97 SPONSORING 7 MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, D.C. 20546-001 


S FUNDING NUMBERS 

C NAS9- 17900 
PR 055-K3-G12 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 

S-671 


10. SPONSORING /MONITORING 
AGENCY REPORT NUMBER 

NASA TP-3207 


14. SUBJECT TERMS 

Strength, Human Modeling, Lean Body Mass, Joint Strength, 
Dynamometer, Force, Torque 

15. NUMBER OF PAGES 
64 

16. PRICE CODE 
A0 4 

17. SECURITY CLASSIFICATION 

18. SECURITY CLASSIFICATION 

19 SECURITY CLASSIFICATION 

20 LIMITATION OF ABSTRACT 

OF REPORT 

OF THIS PAGE 

OF ABSTRACT 


Unclassified 

Unclassified 

Unclassified 



Standard Form 298 (Rev 2-89) 
Prescribed by ANSI Std 239-18 
298-102 


NASA-Langley, 1992 
















